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Abstract 

For a large class of scalar-tensor-like gravity whose action contains nonminimal couplings be¬ 
tween a scalar field (p(x <r ) and generic curvature invariants (R| beyond the Ricci scalar R = R a a , we 
prove the covariant invariance of its field equation and confirm/prove the local energy-momentum con¬ 
servation. These (p{x a ) - % coupling terms break the symmetry of diffeomorphism invariance under a 
particle transformation, which implies that the solutions to the field equation should satisfy the con¬ 
sistency condition R = 0 when <p(x a ) is nondynamical and massless. Following this fact and based on 
the accelerated expansion of the observable Universe, we propose a primary test to check the viability 
of the modified gravity to be an effective dark energy, and a simplest example passing the test is the 
“Weyl/conformal dark energy”. 

PACS numbers: 04.20.Cv , 04.20.Fy , 04.50.Kd 
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1. Introduction 

An important problem in relativistic theories of gravity is the divergence-freeness of the field 
equation and the covariant conservation of the energy-momentum tensor. In general relativity (GR), 
Einstein’s equation G MV = R /n , - \Rf flv = S7tGT i l ' l y ) has a vanishing covariant divergence due to the 
contracted Bianchi identities V /J G /Jv , = 0, which guarantees the local energy-momentum conservation 
WT™ = 0 or s/-gT ( ”^) = 0 for the matter tensor Tf%\ In modified gravities beyond GR and its 
Hilbert-Einstein action, the conservation problem becomes more complicated and has attracted a lot 
of interest. 

In Ref.[l], the generalized Bianchi identities were derived for the Palatini formulation of the non¬ 
linear f(R) gravity, and its local energy-momentum conservation was further confirmed in Ref. [2] by 
the equivalence between Palatini f(R ) and the a> = -3/2 Brans-Dicke gravity. Ref. [3] investigated a 
mixture of f(R) and the generalized Brans-Dicke gravity, and proved the covariant conservation from 
both the metric and the Palatini variational approaches. For Einstein-Cartan gravity which allows for 
spacetime torsion, both the energy-momentum and the angular momentum conservation were studied 
in Ref.[4] by decomposing the Bianchi identities in Riemann-Cartan spacetimes. In Refs. [3, 5, 6, 7], 
the nontrivial divergences were analyzed for the situations where the matter Lagrangian den¬ 

sity is multiplied by different types of curvature invariants in the action. Also, interestingly in Ref. [8], 
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the possible consequences after dropping the energy-momentum conservation in GR, such as the mod¬ 
ified evolution equation for the Hubble parameter, were investigated. 

Besides the covariant invariance V / '7’ / ( ,'''' 1 = 0 for the matter tensor Tffi that has been standardly de¬ 
fined in GR and modified gravities (cf. Eq.(20) below), the conservation problem has also been studied 
for more fundamental definitions of energy-momentum tensors from a wider perspective, i.e. from a 
first-principle approach making use of Noether’s theorem and the classical field theory. For example, 
the Noether-induced canonical energy-momentum conservation for the translational invariance of the 
Lagrangian was studied in Ref. [9] for general spacetimes with torsion and nonmetricity. The conser¬ 
vation equations and the Noether currents for the Poincare-transformation invariance were studied in 
Ref. [10] for the 3+1 and 2+1 dimensional Einstein gravity and the 1 + 1 dimensional string-inspired 
gravity. Also, Refs.[l 1] and [12] extensively discussed the diffeomorphically invariant metric-torsion 
gravity whose action contains first- and second-order derives of the torsion tensor, and derived the full 
set of Klein-Noether differential identities and various types of conserved currents. 

In this paper, our interest is the covariant invariance of such modified gravities whose actions 
involve nonminimal couplings between arbitrary curvature invariants {R.) and a background scalar field 
cf>(x a ). For example, <p(xR) is coupled to the Ricci scalar R = R a a in Brans-Dicke and scalar-tensor 
gravity in the Jordan frame [13], to the Chem-Pontryagin topological density in the Chem-Simons 
modification of GR [14], and to the Gauss-Bonnet invariant Q = R 2 - AR a pR a P + R ( ^p v R' n ‘^ v in the 
Gauss-Bonnet effective dark energy [15]. In theory, one could consider the nonminimal coupling of 
(p(x a ) to an arbitrarily complicated curvature invariant beyond the Ricci scalar. In such situations, 
however, the covariant invariance of the field equation has not been well understood, so we aim to 
carefully look into this problem by this work. Note that it might sound more complete to analyze the 
global conservation V / '(7’ / ( ", !l + t fn ,) = 0 or d M \-^ zr g(T ( ™ > + ? /(V .)| = 0, where t^ v refers to the energy- 
momentum pseudotensor for the gravitational field, but to make this paper more clear and readable, 
we choose to concentrate on the local conservation V /J T ‘"! 1 = 0, while the incorporation of t fJV will be 
discussed separately. 

This paper is organized as follows. In Sec. 2, we introduce the generic class of modified gravity 
with the nonminimal r/TGj-couplings to arbitrary Riemannian invariants [R], calculate the divergence 
for different parts of the total action, prove the covariant invariance of the field equation, and confirm 
the local energy-momentum conservation. Section 3 investigates the reduced situations that the scalar 
field is nondynamical and massless, and derives the consistency constraint R = 0 which suppresses 
the breakdown of diffeomorphism invariance. Finally, applications of the theories in Secs. 2 and 3 
are considered in Sec. 4. Throughout this paper, we adopt the geometric conventions T" y = r 0 ^ , 
R a p yS = dyT*, - d s T^p ■■■ and R )IV = R ,r fJ(rv with the metric signature (-, + + +). 

2. General theory 

2.1. Scalar-tensor-like gravity 

Consider a theory of modified gravity or effective dark energy given by the following action, 

S = J d 4 X \Pg (.%£ + X C + i? NC + Xi) + S m , (1) 

where J2 ?he refers to the customary Hilbert-Einstein Fagrangian density as in GR, 


-5+tE - R , 


( 2 ) 
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while Jz?g denotes the extended dependence on generic curvature invariants R, 


Jz? G = (3) 

Here R = R (g a p , R f , ( ,vfi , ,...) is an arbitrary invariant function of the metric as well as the 

Riemann tensor and its derivatives up to any order. For example, R can come from the fourteen 2 alge¬ 
braically independent real invariants of the Riemann tensor [16] and their combinations, say R a pR a ^ + 
RanpvR a ^ v + Ra/j/ 3 vR a ^R MV , which will yield fourth-order field equations; or differential Riemannian 
invariants that will lead to sixth- or even higher-order field equations, like RV a V ct R+R a ^ v V a VpR^ v . 
In the total Lagrangian density, J2 ?nc represents the nonminimal coupling effects, 

-2kc = ft(0) ■/(*>■■• ,*) , (4) 

where h{<f>) is an arbitrary function of the scalar field (f> = <p{x a ), and /(/?,■■■ , < R) has generic depen¬ 
dence on curvature invariants, with the dots “■ ■ ■ ” in /(/?,•■ ■ , R) and the /(/?,••■ , R) above denoting 
different choices of R. Moreover, the kinetics of (f>(x a ) is governed by 

-2* = -A(<P) ■ V a 0V> - V(0). (5) 

In the (-, + + +) system of conventions, <p(x a ) is canonical if A(cp) > 0, noncanonical if A(cp ) < 0, and 
nondynamical if A{(f>) = 0. 

Finally, as usual, the matter action S m in Eq.(l) is given by the matter Lagrangian density via 

S m - \6nG f d 4 x V 3 # (&V, t/fin, d^m) , (6) 

where the variable t/r„, collectively describes the matter fields, and ij/ m is minimally coupled to the met¬ 
ric tensor g flv . Unlike the usual dependence on d /( (A,,, in its standard form, S£ m = S£ m [g^ v , <A/«, 
does not contain derivatives of the metric tensor - such as Christoffel symbols or curvature invariants, 
in light of the minimal gravity-matter coupling and Einstein's equivalence principle; physically, this 
means .A£ m reduces to the matter Lagrangian density for the flat spacetime in a freely falling local 
reference frame (i.e. a locally geodesic coordinate system). 

To sum up, we are considering the modifications of GR into the total Lagrangian density = R + 
f(R , • • • ,'R)+h((p)-f (R, ■ ■ ■ , < R)-A{(p)-V a (pV a (t)-V((t>)+l67rGJ£ m , which has been rescaled so that the 
numerical coefficient 1 6nG is associated to ,A£ m . It can be regarded as a mixture of the nonlinear higher- 
order gravity «£? = R + f (R, ■ ■ ■ , R) + 1 6nG.2 , lll in the metric formulation for the curvature invariants, 
and the generalized scalar-tensor gravity Jz? = h(cp) •/(/?,••• , < R)~ A(cp) ■ V ( ,</;V"r/>- V(cp) + 1 6nG. r Z J m in 
the Jordan frame for the scalar field. Hereafter we will mainly work with actions, and for simplicity we 
will sometimes adopt the total Lagrangian density in place of the corresponding action in full integral 
form. 


2 When one combines the spacetime geometry with matter fields in the framework of GR. the amount of independent 

algebraic invariants will be extended to sixteen in the presence of electromagnetic or perfect-fluid fields [17]. 
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2.2. Divergence-freeness of gravitational field equation 


2.2.1. Pure curvature parts 

For the Hilbert-Einstein part of the total action, i.e. <She = f d 4 x Jz?he, its variation with 

respect to the inverse metric yields the well-known result &She = f d 4 x sJ—gG^dg^. By the symbol 
= we mean the equality after neglecting all total derivatives in the integrand or equivalently boundary 
terms of the action when integrating by parts, and the Einstein tensor G^ v = - ^Rg fJV respects the 

twice-contracted Bianchi identity = 0. 

For the gravitational action <Sg = f d 4 x x/-g S/fi for the extended dependence on generic Rie- 
mannian invariants, formally we write down the variation as <5 <Sg = f d 4 x yf^g H { ^ } 5g^ v , where H^) 
resembles and generalizes the Einstein tensor by 


r(G) ^ 1 d\f^gf{R,---,R)\ 


(7) 


Due to the coordinate invariance of 5 g, satisfies the generalized contracted Bianchi identities 
[18, 19] 




1 ^ [V=g 

k vFs tig 41v 


( 8 ) 


or just = 0 by the definition of H^\ Similar to the relation G /n . = R IIV - \Rg )JV , one can further 

expand to rewrite Eq.(8) into 


Lv + 2 m MV - ■ • •, m g J = o, 


(9) 


where ()< := df(R, ■ ■ ■ , K)/dR. fin := df(R. ■ ■ ■ ,'R)/dR, and R fn , = (finSR) /Sg MV - note that in the 
calculation of R^ v , fin will serve as a nontrivial coefficient if fin + constant and should be absorbed 
into the variation 6R. when integrated by parts. 


2.2.2. Nonminimal f(x a )-cur\>ature coupling part 

For the componential action <Snc = f d 4 x s[-g ,2nc for the nonminimal coupling effect, formally 
we have the variation &Snc = f d 4 x \/ ir g H ( l ^ C> dg IJV , where 


r (NC) 1 sy=g<f>f{R,--- ,K)\ 

MV sFs 


( 10 ) 


Unlike with Eq.(8), H^ C) does not respect some straightforward generalized Bianchi identities; 
this is because <Snc involves the coupling with the background scalar field (p(x <t ) and is no longer 
purely tensorial gravity. Thus, we will analyze the divergence of H^ C) by the diffeomorphism of 
<5nc- 

Consider an arbitrary infinitesimal coordinate transformation x 1 ' i-» x IJ + 6x IJ , where 8xP - k fl is 
an infinitesimal vector field that vanishes on the boundary, = 0 | an , so that the spacetime manifold 
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is mapped onto itself. <Snc responds to this transformation by 


&S N c = - J' d 4 x h(<p) ■ df, [kf V 3 ?/ (R, ■ ■ ■ , 7?)] 

= J d 4 x V=^ f(R, ■ • ■ , K) ■ (hfdrf) kP, 


( 11 ) 

( 12 ) 


where h,/, := dh((p)/d<p. For Eq.(ll), one should note that 4>(x <r ) acts as a fixed background, as it 
only relies on the coordinates (i.e. spatial location and time) and is independent of the spacetime 
metric; moreover, the coordinate shift x /J i-» W + kP is a particle/active transformation, under which 
the dynamical tensor field g^ v and thus x[-g f (/?, ■ •■ , R) transform, while the background field tp{x a ) 
and the coordinate system parameterizing the spacetime remain unaffected [20], 

Under the particle transformation x /J x 1 ' + k>\ the metric tensor varies by g^ v i-» g /JV + 6g M y with 
5g^ v - £%gnv - -V ' M k v - Vyk^, and therefore g^ v i-> g MV + dg^' with 5g^ v = -£^ v = + V v k p . 

Recalling the definition of H ( ^ C) with H ( ^ C) being symmetric for the index switch p <h> y, one has 


&Snc = 2 J d 4 x V~g H i f ^ C) V tJ k v = -2 r d 4 x V=^ (v^ C) )k v , (13) 

Comparing Eq.(12) with Eq.(13), we conclude that H^ c> has a nontrivial divergence for (p(x rv ) + 
constant, and 

V^ C) = /(/?,••• ,K)-htdy4>. (14) 

In fact, Eq.(14) reflects the breakdown of diffeomorphism invariance in the presence of a fixed back¬ 
ground scalar field. As a comparison, it is worthwhile to mention that under an observer/passive trans¬ 
formation where the observer or equivalently the coordinate system transforms, both the tensor fields 
and the background scalar field will be left unchanged, so the symmetry of observer-transformation 
invariance continues to hold [20]. 


2.2.3. Purely scalar-field part 

Next, for the purely scalar-field part S ( i, = f d 4 x sf^g with the variation dS ( t, = f d 4 x rigdg^ v , 
explicit calculations find 

rig = ~m ■ V„0V v 0 + ■ V ff 0V"0 + V(0)) g^v ■ (!5) 

Taking its contravariant derivative, we immediately obtain the nontrivial divergence 

V*rig = ^ 0V ff 0 + 24(0) • Q0 - V (! ) • V v 0, (16) 

where A,/, := dA((p)/d(f>, V ( i, := dV((p)/d(f>, and □ denotes the covariant d’Alembertian with D0 = 
g a/3 V a Vp<P = -^= da{-\f~g g a ^dpf). On the other hand, extremizing the entire action Eq.(l) with respect 
to the scalar field, i.e. SS/6<f> = 0, one could obtain the kinematical wave equation 


24(0) ■ D0 = -/ (R, ■ • • , K) ■ h+ - 4 0 • V a 0V“0 + . (17) 

We regard it as “kinematical” because it does not explicitly relate the propagation of <p(x a ) to T (m> - 
g^T^ for the matter distribution, while the “dynamical” wave equation can be obtained after comb- 
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ing Eq.(17) with the trace of the gravitational field equation. Substitute Eq.(17) into the right hand 
side of Eq.(16), and it follows that 

( 18 ) 

which exactly cancels out the divergence of H { ^ c> in Eq.(14) for the nonminimal-coupling part 5 nc- 
2.2.4. Covariant invariance of field equation and local energy-momentum conservation 

To sum up, for the modified gravity or effective dark energy given by Eq.(l), its field equation 
reads 


C* 4- Z/(^) _i_ Z/(NC) rr(0) _ Q 7r /^T'( m ) 


(19) 


where, unlike G jiv and H^J. the exact forms of {will not be determined until the con¬ 
crete expressions of {2£g, 2z?nc} or {/(/?, • • • ,7?) ,f(R, ■ ■ ■ ,7?)} are set up. In Eq.(19), the energy- 
momentum tensor T { ™ ] is defined as in GR via [21] 

dS,n = ~ X 16 nG J d 4 x T^6g MV with T ( ™ } := —- , ( 20 ) 


with S m rescaled by I 6ttG in Eq.(6). Instead of the variational definition Eq.(20), it had been suggested 

m =0 for the i p m 


that could be derived solely from the equations of motion „ . - V 


L U v 


dfi„ 


(Xfiitfim) 


field in [g MV , ip m , [22]; however, further analyses have shown that this method does not hold 

a general validity, and Eq.(20) remains as the most reliable approach to 1 [23]. 

Adding up the (generalized) contracted Bianchi identities V^G^y = 0 and Eq.(8), and the nontrivial 
divergences Eqs.(14) and (18), eventually we conclude that the left hand side of the field equation (19) 
is divergence free, the local energy-momentum conservation = 0 holds, and the tensorial 

equations of motion for test particles remain the same as in GR. 

In fact, the matter Lagrangian density f£ m = S£ m [g^v, ip m , d^ipm) is a scalar invariant that respects 
the diffeomorphism invariance under the particle transformation x IJ h-> x IJ + k 1 ', and Noether’s conser¬ 
vation law directly yields 




xFg 


( 21 ) 


which can be recast into = 0. That is to say, under minimal geometry-matter coupling with 

an isolated d£ m in the total Lagrangian density, the matter tensor in Eq.(20) has been defined in 
a practical way so that T ( ™ ] is automatically symmetric, Noether compatible, and covariant invariant, 
which naturally guarantees the local conservation = 0. In this sense, one can regard the 

vanishing divergence W (G flv + + H ( ^ C> + = 0 for Eq.(19) to either imply or confirm the 

conservation WTffi = 0. 

One should be aware that in the presence of nonminimal gravity-matter couplings, like 7? • 
or more generally F (R,- ■ ■ , K) ■ F£ m in the total Lagrangian density, the divergence V /J T' t 'f ] becomes 
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nonzero as well and obeys the relation = F (R,- ■■ , 'R) 1 • (FF m g MV - Tffi'j • ¥ IJ F (/?,••• , R) 

instead [5, 6, 7], which recovers the local conservation V / '7’ / ( ,'" ) = 0 for F (R, ■ ■ ■ ,'R) = constant. 

Also, at a more fundamental level, the in Eq.(20) for GR and modified gravities, though 
practical with all desired properties, is not defined from the first-principle approach, i.e. directly 
from symmetry and Noether’s theorem in the classical field theory. In this larger framework, the 
in Eq.(20) is often referred to as the Hilbert energy-momentum tensor: it symmetrizes the canonical 
energy-momentum tensor of translational invariance by adding a superpotential term, and it is a special 
case of the Belinfante energy-momentum tensor that minimally couples to gravity [24]. 

3. Nondynamical massless scalar field 

3.1. Nondynamical massive scalar field 

Due to the T(0)-dependence in S ( t„ its Lagrangian density becomes 2z?^ = —V((p) when AUp) = 0; 
considering that V{(p) is usually related to the mass of the scalar field in cosmology and high energy 
physics, we will call cp{x a ) nondynamical and massive for the situation A(cp) = 0 and VUp) + 0. 
As such, instead of producing a prorogation equation ucp. the extremization SS/5(p = 0 leads to the 
following constraint for the potential V{(p)\ 




( 22 ) 


In the meantime, Eqs.(15), (16), and (17) reduce to become 

= \ V (0) 8rv and = \ V ( ,NA = \j(R, ■,‘R)- h 4 N A . (23) 

Thus, for a nondynamical yet massive scalar field, WHffi can still balance the nontrivial divergence 
V^ C) of the nonminimal 0(x Q ')-eurvature coupling effect, while the potential or the mass of the 
scalar field is restricted by the condition Eq.(22). 

3.2. Nondynamical massless scalar field 

Within the situation A{(p) = 0, it becomes even more interesting when the potential vanishes as 
well in Eqs.(5), (15), (16), and (17); we will call the scalar field nondynamical and massless 3 for 
= 0 = V(cp). With = 0, the total action simplifies into 

S ~ J' d 3 4 x s[-g(R + .-^G + -^nc + 1 . (24) 

Since H^J = 0 and V = 0, the divergence V 1 'H 1 ^ 1 for the nonminimal coupling part as in 
Eq.(14) can no longer be neutralized. Instead, with = 0, the generalized contracted Bianchi 

identities Eq.(8), and the covariant conservation V / '7’ / ! "' , = 0 under minimal geometry-mater coupling, 
the contravariant derivative of the field equation G t , v + Flffi + H ( ^, C) = 8nGT^ forces to 


3 We simply use “massive” and “massless” to distinguish the situation V(p) t- 0 from V(cf>) = 0 when the scalar field is 

nondynamical. However, we do not follow this usage to call <p(x a ) “dynamical and massless” when {A{(p) + 0, V(cp) = 0), as 
it sounds inappropriate to from the spirit of relativity. 
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vanish. Together with Eq.(14), this implies that to be a solution to the gravity of Eq.(24), the metric 
tensor g fJV must satisfy the constraint 

/(/?,•■■ , R) = 0 for (f){x a ) + constant. (25) 

Since the nonzero divergence V 1 ' = ~^f(R, ■ ■ ■ ,'R) ■ h,/,d v (l) measures the failure of diffeomor- 
phism invariance in the componential action 5 no the consistency condition Eq.(25) indicates that the 
symmetry breaking of diffeomorphism invariance is suppressed in gravitational dynamics of Eq.(24). 

Here one should note that the variation 6S/6(p = 0 yields the condition f (R, - ,‘R) ■ ft# = 0, 

which also leads to f (R, ■■■ ,‘R) = 0 if the scalar field is nonconstant. In addition, the constraint 
/(/?,•• ■ ,R) = 0 does not mean H^ C) = 0 or the removal of 2 ??nc from the action Eq.(24). This can 
be seen by an analogous situation in GR: all vacuum solutions of Einstein’s equation have to satisfy 
the condition R = 0, but the GR action S = f d 4 x -\f~g (R + 1 67 rG 2 z? m ) still holds in its standard form. 

After 2 z?g and .iAic get specified in Eq.(24), how can we know whether it yields a viable theory or 
not? In accordance with Eq.(25), we adopt the following basic assessment. 

Primary test: For the action Eq.(24) to be a viable modified gravity or effective dark energy carrying a 
nondynamical and massless scalar field, an elementary requirement is that the function f (R, - ■ ■ , R) in 
J^nc vanishes identically for the flat and accelerating Friedmann-Robertson-Walker (FRW) Universe 
with the metric 

3 2 

ds 2 = -dt 2 + a(t) 2 ^ (d x ') and d(t) > 0, (26) 

1=1 

where a{t) is the cosmic scale factor and the overdot means derivative with respect to the comoving 
time. 

This primary test is inspired by the fact that the observable Universe is homogeneous and isotropic 
at the largest cosmological scale, and the discovery that the Universe is nearly perfectly flat and cur¬ 
rently undergoing accelerated spatial expansion. These features have been extensively examined and 
received strong support from the surveys on the large scale structures, the expansion history, and the 
structure-growth rate of the Universe, such as the measurements of the distance modulus of Type la 
supernovae, peaks of the baryon acoustic oscillation, and temperature polarizations of the cosmic mi¬ 
crowave background. Clearly, the primary test is updatable and subject to the progress in observational 
cosmology. 

3.3. Weyl dark energy 

Following the primary test above, one can start to explore possible modifications of GR into the 
total Lagrangian density 2z? = R + f(R. ■ ■ • ,'R) + h(<p)f(R, • • • , R) + 16zr GJ£ m and then check the 
consistency condition /(/?, ■ • • ,R) = 0 under the flat FRW metric Eq.(26). In the integrand of the 
Hilbert-Einstein action for GR, the Ricci scalar R is the simplest curvature invariant formed by second- 
order derivatives of the metric; similarly, we can start with the simplest situation that /(/?,■■■ ,'R) is 
some quadratic Riemannian scalar. One possible example is the square of the conformal Weyl tensor 
Caf}y6 — Rafiyd 3" 2 — SayRflS 3" Spy^ad ~ Spd^ay^ 3" g {SaySpS ~ SaSSpy) which is the totally 

traceless part in the Ricci decomposition of the Riemann tensor. In this case, we consider the action 

S C 2 = f d 4 x (R + '!<]> C 1 + 1 6nGXn) , 


(27) 


where y ± 0 is a coupling constant, and 


/~>2 ._ /^auBv — ^ r>2 ^ r> r>wv i d p/iQ'v/? 

^ •— '-'anfiv'-' — 'LKiavK “l" KuavBK 






(28) 


It is straightforward to verify that C 2 = 0 for arbitrary forms of the scale factor a{t ) in the flat FRW 
metric. We would like to dub Eq.(27) as the “Weyl dark energy” or “conformal dark energy”. The 
field equation is R jn , - \Rg IJV + y H^, 1 = StiGT^, where 

Hfiv = ~ ~ 2 R^ Rav RfiaftyRy " /J> j + | (&«,□ - V„V v ) (<f>R) - 2n (cpR^) 

+ 2V„V v (0*/) + 2(« pR v a ) - 2g pv V £1 ,V p (<f>R°#) + 4V^V“ (cpR a ^ v ) , (29) 

and according to Eq.(14), its covariant divergence is 

(30) 


v^ 2) = -|c 2 v v 0. 


The Weyl dark energy S C 2 , where the scalar field is nondynamical and massless, can be generalized 
into the dynamical case 


S = J d 4 x V=g (r + h(<p)C 2 - 4(0) ■ y a (p^ a (p - V(<p) + 16 jtG ^, (31) 

for which the constraint C 2 = 0 is no longer necessary and should be removed. 

The complete validity of the Weyl dark energy S c 2 or its extension Eq.(31), including the value 
of the coupling strength y in •Sc 2 ’ should be carefully constrained by the observational data from 
astronomical surveys. Following the field equation of S C 2, consider a C 2 CDM model (i.e. C 2 cold 
dark matter) for the Universe instead of ACDM. Then the first Friedmann equation under the flat FRW 
metric reads 




( a o \ 4 

+ PK>(-J + P C 2 > 


(32) 


where the densities of nonrelativistic matter pm( 0 and relativistic matter p r (t) have been related to their 
present-day values pmo and p,-o via by the continuity equation p + 3 Hp{ 1 +w) = 0, with the equation of 
state parameters being wm = 0 and w r = 1/3, respectively. Also, H := d/a is the evolutionary Hubble 
parameter, and p c 2 denotes the effective energy density of the Weyl dark energy, 


. a a . d 2 a . a J . d J a a a J 

Pc 2 = Y ( 5 0- 2 0— - - 0— - 20— + 5(f>— - + 2 0— 

a a a 1 a a* a 3 a a a* 


. a 3 


. d 5 


. 2 

a a u u a cr cr uu a i 

■40- — - 40-+ 40— — - 60—- + 80— . 

a n L n n n n L n L n rn J 


a a z 


a a a 

) - 

a a a 


2 • 2 
a a 


cr a ^ 


a a 
b —- 
a a 4 


a 


(33) 


Employing the cosmological redshift z '■= oo/a - 1 as well as the replacements a/a = H + H 2 and 
a /a = H + 3 HH + H \ Eq.(32) can be parameterized into 


H(z-,Ho,p) = H 0 


^M0( 1 


+ z) 3 + fl r o(l + z) 4 + U ( 


te 2 


(34) 
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where Hq represents the Hubble constant H(z = 0), £2mg = SkGpmo/OH 2 ), £2,-o = SnGp r o/(3H^), and 

n C 2 = ( H + H 2 ) - 2 [h + H 2 f -H 2 -2H 4 \ + (f>H (5-4 H- AH 2 ) ( H + 3 HH + H 3 ) 

+ 8 (pH 6 + (p[H + H 2 )\ (2 H + 2H 2 + AH 2 ) [ti + H 2 ) - AH 2 - 6H 4 ] J . (35) 

Typically, we can use the Markov-Chain Monte-Carlo engine CosmoMC [25] to explore the parameter 
space p = HI mo, £2,o, y) for the Weyl dark energy S c i, and find out how well it matches the various 
sets of observational data. This goes beyond the scope of this paper and will be analyzed separately. 

4. Applications 

4.1. Chern-Simons gravity 

The four-dimensional Chem-Simons modification of GR was proposed by the action [14] (note 
that not to confuse with the traditional gauge gravity carrying a three-dimensional Chem-Simons term 
[26]) 


S C s = J" d 4 x V 3 ! + Y0^== + 167rG ^„j . (36) 

The scalar field (p = (p(x n ) is nonminimally coupled to the Chern-Pontryagin density *RR := * R a / 3 -ysR a ^ y6 
^aPnvR^ yS R aPy6 , where *R a p yS := ^e a p MV R MV y5 is the left dual of the Riemann tensor, and e a p MV rep¬ 
resents the totally antisymmetric Levi-Civita pseudotensor with eoi23 = sf~g and e 0123 = 1/ xf-g. 
The field equation reads R jiv - \Rgfj V + yH { ^ = 8nGT ( ™\ where // / ( , ( 1 ; P| = —== 6{ ^ R , R} collects the 
contributions from the ^(.r Q ')-coupled Chem-Pontryagin density, 

V=2 < p) = 2 • (v/jVV v + + 2 dadrf ■ (*R a /y + *R a / fJ ) . (37) 

According to the general theory in Secs. 2.2 and 3.2, the Chern-Simons gravity Eq.(36) involves 
a nondynamical and massless scalar field. Identifying /(/?,-■• ,R) as RR/ H~g and with h((p) = yep 
in Eqs.(14) and (25), we obtain the divergence 


V^ P) = -l^.d v cp, (38) 

2 V=g 

as well as the constraint *RR = 0 for nontrivial (p(x a ). It can be easily verified that RR vanishes for the 
flat and accelerating FRW Universe, and thus passes the primary test in Sec. 3.2. Also the condition 
*RR = 0 only applies to the action Eq.(36), and is invalid for the massive Chern-Simons gravity 
j£f = R + Y0~|== _ + 167r GJ£ m or the dynamical case 2zf = R + y<P~j= ~ A((p) ■ V a <pV a <p + \6nGH^ m . 

4.2. Reduced Gauss-Bonnet dark energy 

The Gauss-Bonnet dark energy was introduced by the action = J d A xx[~g{R + h{ep)Q - 

(t (p¥ ir (p - V(<p) + 16/rG2^„j [15], where A £ {±1,0}, and the scalar field is nonminimally coupled to 
the Gauss-Bonnet invariant Q := ( \ e af3y[,R y ^) ■ ( I e, l ( p(T R fim, d) = R 2 - 4R a pR a P + R a ^ v R a ^ v ■ If <P(x ,r ) 
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is nondynamical with X = 0, the action .S^ reduces to become >S[?g = f d 4 x \R + h{<p)Q - V(cp) + 

16 7rGJZ? m 'j, and according to Eq.(22) with /(/?,•■• /R) identified as the Gauss-Bonnet invariant, V(cp) 

has to satisfy the constraint V,/, - Q h r/) . Moreover, the nonminimally coupled hi<p)Q part in <S[ig and 
«S[?g contributes to the field equation by 

//^ B) = 2 R ( gfIV D - V^V v ) h + 4 R/VJyh + 4R v a V a V M h 

- 4R MV ah - 4gfj V ■ R^VaVph + 4R a ^V a h , (39) 

where, compared with the original field equation in Ref.[15], we have removed the algebraic terms in 
H ( fi B) by the Bach-Lanczos identity \Q gflv = 2RR^ - 4R^R av - 4R cmPv R a P + 2 R lufiy R v apy [27]. The 
divergence of H ( ^ B \ in accordance with Eq.(14), reads 

v^° B) = - l -g.h' P d v( p. (40) 

However, it would be problematic if one further reduces «S[fg into 

^gb = f d4x [ R + K4>)G + 167tGJ^„) , (41) 

where cp{x a ) is both nondynamical and massless. The metric tensor has to satisfy Q = 0 to be a solution 
to the field equation R /n , - \R gflv + H^ B> = 8 nGT^ for the reduced Gauss-Bonnet dark energy S £?g. 
For the flat FRW Universe with the metric Eq.(26), the Gauss-Bonnet invariant is 

G = 24^ , (42) 

a i 

and thus Q vanishes only if the Universe were of static state (a = 0) or constant acceleration (a = 0). 
Hence, the constraint Q = 0 for S G g is inconsistent with the cosmic acceleration, which indicates that 
unlike S j.lg and S G g, S^g is oversimplified and can not be a viable candidate of effective dark energy. 

4.3. Generalized scalar-tensor theory 


Since 5 he and Sq in Eq.(l) respect the diffeomorphism invariance and the (generalized) con¬ 
tracted Bianchi identities, in this subsection we will ignore them and focus on the following scalar- 
tensor-type gravity in the Jordan frame: 


<Sst = 



V=g(/(R, <P) + -S?nc + £’</>+ 167rG«Sf m j 

V=g (f(R, <f>) + h(<P) .f(R,...,<R)- A((p) • V a cPV a <P - V(cp) + 16 nG^n), 


(43) 


where f(R, <P) is a hybrid function of the Ricci scalar and the scalar field. f(R. (p) contributes to the 
field equation by 


H, 




VV 


1 d{s/^gf(R,(p)) 

xFs 


~ X.f( R > <P ) ' g[iv + fRR,v + - V^V v ) f R . 


(44) 
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where J'r = Jr(R, (f>) = df(R,(p)/dR. With the Bianchi identity V /J - ^Rg^ = 0 and the third- 
order-derivative commutator (V v n - aV v ) fa = explicit calculations yield 

= _l /<4 . Vy0) (45) 

where = f/,(R, (f>) = df(R,cf>)/d(t>. On the other hand, the kinematical wave equation 6S$j/d(p = 
0 reads 24(0) • n<p = -fy - f ■ h# - A<p ■ V a (f>V a (p + V^, which recasts the divergence = 

• V a cl>V a (t> + 24(0) ■ D0 - V$) ■ V v 0 as in Eq.(16) into 

= \[u + f(R, ■■■,*)■ h^j ■ V v 0 . (46) 

Hence, with Eqs.(14), (45) and(46), we immediately learn that the field equation H + H^ C) + 
H ( J’y = SnGT { "y for the scalar-tensor-type gravity «Sst is divergence free. By the total Lagrangian 
density for the sake of simplicity, the concretization of Eq.(43) includes, for example, standard Brans- 
Dicke gravity Jz? = (f>R - V a cf>V a (p + I67rj£ m (where Newton’s constant G is encoded into 0 1 in 

the spirit of Mach’s principle) [13], generalized Brans-Dicke gravity J§? = <pR - V a (f>V a (p - V((p) + 

l6nGJt? m with a self-interaction potential, Lovelock-Brans-Dicke gravity S£ = 0 |a’ + ~^=*RR + b(j j- 

^V a 0V“0 - V((l>) + 1 6nG.£m [27], Lovelock-scalar-tensor gravity 2zf = f\((f>)R + yj=RR + ,h(<j>)Q - 
a>(0) • V ff 0V"0 - V(0) + 167rG2z^„ [27], minimal dilatonic gravity Jz? = <pR - 2AG(0) [28], Gauss- 
Bonnet dilatonic gravity 2zf = R - V a (/>V a (p + e~ y ^Q or Jz? = e~ r ^(R - V a 4>V a (t> + Q) motivated by the 
low-energy heterotic string theory [29], and the standard scalar-tensor gravity Jz? = F(cf>)R - Z( 0) ■ 
V a 0V“ - V(0) + 1 6nGJ£ m [30]; all these examples satisfy the local energy-momentum conservation 
= 0 and have divergence-free field equations. 

4.4. Hybrid metric-Palatini f(R) gravity 

So far we have been using the metric formulation for the curvature invariants; however, the local 
conservation can be proved for the Palatini or hybrid metric-Palatini f(R) gravity without referring to 
the Palatini formulation of the (generalized) Bianchi identities. Consider the following hybrid metric- 
Palatini f(R) action 


S Hf = f dA x^g(R + f(R) + 167rGJz? m ), (47) 

A A A /V, ZV 

where R is the usual Ricci scalar for the metric g t , v , while R = R(g, T) = g llv R llv (T) denotes the Palatini 
Ricci scalar, with the Palatini Ricci tensor given by R^ v (t) = R a Mav (T) = d a t - d v t^ + f“ f^ v - 
Variation of <Sh/ with respect to the independent connection f“ v yields V„( yj—gffcg^) = 0, 
where V is the covariant derivative of the connection and fj { := df(R)/dR. Thus, V is compatible with 
the auxiliary metric f fi g IJV =: g MV , as \p£g IJV = xFgffrg^. Relating g fJV to g^ v by the conformal 
transformation g flv g )Jv . and accordingly rewriting R fiv and R in the metric formulation, one could 
find that .Sj] j is equivalent to [31] 

= f d 4 x ^>(R + $R+ V a <pV a <f> - V(0) + 16;rGJz^), (48) 
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where 0(x“) = fj t (R) and U(0) = f^R-f(R). is just the mixture of GR and the o»bd = -3/2 Brans- 
Dicke gravity. Recall that Eq.(43) has employed the generic function f(R. 0) for 5 st> which includes 
the hybrid situations like /(R, 0) = R + (pR. Hence, following Sec. 4.3, it is clear that the hybrid 
scalar-tensor gravity and thus the hybrid metric-Palatini f(R) gravity have divergence-free 
field equations and respect the local energy-mentum conservation. 


5. Conclusions 


In this paper, we have investigated the covariant invariance of the field equation for a large class of 
hybrid modified gravity 55 - R+f (R, ■ ■ ■ ,R) + h(<p) • / (R, • • ■ /R) - A((p) ■ V a (pV cr <p-V(tf>) + 16 7tG55 m . 
For the four components «£he = R , 2z?g = / (R, • • • , R), 56 nc = K<P) ■ / (R, • ■ • ,R), and 55$ = -4(0) • 
V„0V"'0 - V((p), we have calculated their contributions {to the gravitational 
field equation along with the respective divergences, which proves the divergence-freeness of the field 
equation (19) and confirms/proves the local energy-momentum conservation under minimal gravity- 
matter coupling. 

H^ C) and H ( J’J fail to obey the generalized contracted Bianchi identities due to the presence of 
the background scalar field 0(jt“), but fortunately, the two nontrivial divergences and V /J H^’J 

exactly cancel out each other. When cp{x a ) is nondynamical and massless, i.e. 4(0) = 0 = V(0), 
the divergence V 1 'H l ^ C> = -\f (R, ■ ■ ■ ,R) ■ h,! t d v (p is forced to vanish, which implies the constraint 
/(/?,•■• . R) = 0 for nonconstant (p{x a ). We have suggested a primary viability test for the gravity 
55 = R+f (R, - , ( R)+h{(f>)-f (R, - , R)+ l6nG55 m by requiring that /(/?,■•■ , R) vanishes identically 

for the flat and accelerating FRW Universe, and a simplest example is the Weyl dark energy 55 = 
R + ycpC 2 + l6nG55 m . 

With the general theory developed in Secs. 2.2 and 3.2, we have considered the applications to the 
Chern-Simons gravity, Gauss-Bonnet dark energy, and various (generalized) scalar-tensor gravities. 
In fact, the theory 2 z?st = f(R, 0) + h(<p) •/(/?,••• , R) - 4(0) • V (! ,0V Q '0 - V(0) + l6nGJ5 m in Sec. 4.3 
can be further extended into ,£?est =/(/?,••• , R, 0) - 4(0) • V„0V"0 - VUb) + \6nG55 m , for which 
we conjecture that the covariant conservation still holds, with 


H 


(/) 


11V 


1 ,K0)] 


and WH$ = --MR, ■ • ■ , <R, 0) ■ V v 0. 


(49) 


However, this divergence relation has not yet been proved in this paper, and we hope it could be solved 
in future. 

In prospective studies, we will take into account the existent candidates of the energy-momentum 
pseudotensor t MV for the gravitation field (cf. Ref. [32] for a review), and discuss the global conservation 
W + (m,) = 0. Also, we will make use of more fundamental definitions of the energy-momentum 
tensor, and look deeper into the conservation problem in modified gravities from the perspective of 
Noether’s theorem and the classical field theory. 
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